Exact AdS/CFT spectrum: 
Konishi dimension at any coupling 
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We compute the full dimension of Konishi operator in planar N=4 Super Yang-Mills (SYM) theory 
for a wide range of couplings, from weak to strong coupling regime, and predict the subleading terms 
in its strong coupling asymptotics. For this purpose we solve numerically the integral form of the 
AdS/CFT Y-system equations for the exact energies of excited states proposed by us and A. Kozak. 
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INTRODUCTION 

Four-dimensional Yang-Mills theories are at the heart 
of modern high energy physics, describing all fundamen- 
tal interactions except gravity. Nevertheless, in spite of 
considerable efforts during almost 40 years, we still don't 
have a satisfactory quantitative description of the most 
interesting YM theories, such as QCD, in the region of 
intermediate and strong couplings. The low energy quan- 
tum dynamics of YM field is mostly known only from 
computer simulations of lattice YM theories. A few im- 
portant exact results concerned the topological, BPS sec- 
tors of N=l,2 SYM were obtained. 

Recently, when the hopes on complete exact 4D solu- 
tions in particular, for the quantities given by nontrivial 
4D Feynman series seemed to start waning, N=:4 super- 
symmetric Yang-Mills theory gave us serious hopes for 
a better understanding of the dynamics of strongly in- 
teracting 4D gauge theories. Due to the AdS/CFT cor- 
respondence [1], as well as to the quantum integrability 
discovered on both sides of this duality in the planar limit 
(when the number of colors N ^ oo with the 't Hooft 
coupling A ~ gy-f^N fixed) [2-11], we acquire little by 
little tools for the study of the most important quanti- 
ties in N=4 SYM, such as the spectrum of dimensions 
A(A) of local operators as functions of A - the scale in- 
dependent coupling constant in this superconformal 4D 
theory. The weak coupling behaviour (A — )■ 0) of A(A) is 
given by Feynman perturbation theory. The dual string 
worldsheet cr-model with the coupling g = •\/A/47r allows 
to find the strong coupling asymptotics of various dimen- 
sions. Integrability allows us to connect the two regimes. 
In particular, the asymptotic Bethe ansatz (ABA) of [13] 
gives us the asymptotic spectrum of single trace operators 
containing an asymptotically large number of elementary 
fields. 

However, for short operators, such as Konishi operator 
Tt[D,ZY [12] [39], the calculation of anomalous dimen- 
sions is still an interesting and important challenge. 
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Figure 1: Numerical solution of exact finite size integral Y- 
system equations for the Konishi dimension Aj<-(A) in a wide 
range of 't Hooft couplings A, compared to the asymptotic 
Bethe ansatz curve and to the predicted large A asymptotics 
Ak(A) ~ + IjX^I^ obtained by fit. 



Recently we proposed the Y-system for the planar 
AdS/CFT [14], a set of functional equations defining the 
anomalous dimensions of all operators of planar N=4 
SYM theory at any coupling. The integral form of the 
y-system for excited states in SL{2) sector, including 
the one corresponding to Konishi operator, was presented 
in [15]. The integral equations for the BPS vacuum en- 
ergy were independently obtained in [16, 17] by the ther- 
modynamical Bethe ansatz (TBA) technique based on 
the dynamics of bound states [18-24] (see also [9, 25]) 
of the mirror theory [26, 27]. The solutions of the in- 
tegral Y-system are also solutions of the functional Y- 
system [15-17, 28]. The combination of functional and 
integral versions of the Y-system appears to be quite ef- 
ficient to compute numerically the exact spectrum. In 
this work, we use the functional form of the Y-system to 
derive the large volume (L) [40] asymptotic solution and 
then, departing from it, we solve the integral form of the 
Y-system iteratively. As a demonstration of the power of 
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Konishi state 




Figure 2: Plot of Ak{X) — 2A^^'' from the numerical data 
compared with the Bethe ansatz prediction and some fits. The 
fits in this plot are done assuming the asymptotics (A) = 
2A^/'' + 2/A^/'' + .... 



our method, we calculate numerically the dimension of 
Konishi operator in a wide range of the 't Hooft coupling 
covering both the weak and strong coupling regimes. The 
results appear to be quite satisfactory: we manage to 
compute the dimension of Konishi operator in the inter- 
val of couplings < A < 700 and to confirm, within the 
precision of our numerics, all the existing data concern- 
ing this quantity: the perturbative results [41] up to 4 
loops (up to A''^ terms) [29-32] and the large A asymp- 
totics 2A^/'' matches the prediction of [33] for the lowest 
level of the string spectrum. Fitting our numerical data 
with A > 60 we find (with uncertainty in last digit) 



Ak = 2Ai/** ( 1.0002 
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(1) 

The leading term reproduces indeed the expected large 
A behavior within our numerical precision. It was also 
argued in [34] that the subleading coefficient ought to 
be integer (the next corrections could be transcendental 
[35]). Indeed, our numerics seems to indicate that Ak = 
2A^/'' + ^j}}!^ + . . . thus predicting the value of this 
integer! [42] We also obtained predictions for two further 
subleading corrections (with a lower precision of course) . 
Our results are represented in Figs.l and 2. 



Y-SYSTEM FUNCTIONAL AND INTEGRAL 
EQUATIONS FOR ADS/CFT 

The Y-system defining the spectrum of all local oper- 
ators in planar AdS/CFT correspondence reads [14] 



^a,s^a,s _ (1 + ^a,s-|-l)(l + 5^a,s-l) 
i^a+l.,s5^a-l,s (1 + ^a+l,s)(l + ^a-l^s) 



(2) 



where the functions Yq ,j(w) correspond only to the nodes 
marked by O, ®, ®, Q, • on Fig. 3 (we will use these nota- 



- -O: 



nr-4 4 4 

A— i 

-3 3 3 

-2 2 2 



-a 



s 
> 



Figure 3: T-shaped domain (T-hook) [36]. It defines the in- 
teractions between Y's in the AdS/CFT F-system equations. 



tions for F- functions in what follows). The one particle 
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energy at infinite length e(°^(u) 
fined through the Zhukovski map x{u) + 

= f{u ± m/2), /± = f{u ± i/2) for any function 
f{u) . A solution of K-system with a given set of quantum 
numbers defines the energy of a state (or dimension of an 
operator in N=4 SYM) through the formula (3) where 
the Bcthe roots uj arc fixed by the exact Bethe ansatz 
equations Ym-^{uj) = — 1. In this paper we restrict our- 
selves to the integral form of the F-system for the SL{2)- 
excited states obtained in [15]. Furthermore we focus 
ourself for simplicity on the Konishi operator where we 
have only two roots ui = —U2 which we can encode into 
the "Baxter functions" R^^\u) = {x{u) — xf){x{u) — x^) 
and their complex conjugates B{u) = R{u) where xf2 = 
x{ui^2 T i/2) with the physical choice of branches, such 
that ]a;(w)] > 1, while for the free variable x(u) we 
should use the mirror kinematics which corresponds to 
the branches where Im {x{u)) > [27]. Unless it is ex- 
plicitly said otherwise, we will be always using this latter 
choice in what follows. With the mirror choice, x{u) has 
a semi-infinite cut for u E (— oo, — 2g) U {2g, +oo). The 
energy of the Konishi state is computed from 

oo poo 1 

AK = 2 + 2e(i)(«i) + ^ / -^a„e('')log(l + y,J ,(3) 
where the integral equations defining Y^^ read [15] 

log Y^ = Kra-l * log(l + l/ro,J/(l + Y[^J 



+ 7e("'") *log(l + y.,J + log■ 
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log Yq = Mnrn * log(l + Y,„^ ) - Kn-l^\og{l + Y^ 



log(l + Yq^ 
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logi^o„ = /•i:„-i.,„-^i*log(l + l/YoJ+-ft:„_i®log(l + Yg,) 
log Y.„ = r„™ * log(l + Y.„) + 27e("°)®log(l + Y^) 

+ 2Nnm * log(l + Yq ) + . (4) 
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We use here the kernels and sources defined in [15] and 
presented in the appendix for completeness. The inte- 
grals in convolutions K * f ^ J dvK(u,v)f(v) go along 
the real axis, but slightly below the poles in the terms 
involving Yq^ (due to the last term in the corresponding 
integral equation, see [15]). The convolutions ® should 
be understood in the sense of a B-cycle (see [15]), e.g. 
7^("0)®log(l + Y^) stands for 



2g 



dv 



2s 



7^("") log(l + Y^) - log(l + l/Fe) 



where ^ is the analytical continuation of Y^ across the 

cut u € (— oo, — 2(7) U {2g,+oo). Summation over the 
repeated index m is assumed with m > 2 for £l±m and 
0-|-,„, and m > 1 for •„!. 

A remarkable feature of all these equations, crucial for 
the success of our numerics and noticed in [15], is the 
reality of all y-functions in the integral equations. 

EXACT BETHE EQUATIONS 

The Y-system integral equations for the functions y^.s 
need to be supplemented by the exact Bethe equations 
Y^^{uj) = —1 which reproduce the asymptotic Bethe 
equations of Beisert and Staudacher in the large L limit 
[14]. To use this equation, we need to analytically con- 
tinue the last of eq.(4) in the free variable u from mirror 
to physical plane and then evaluate it at it = mi. After 
some manipulations with contours wc find 

iogi^,pf(^i)=ri™*iog(i+y,„)+iogZQ>i) (5) 

+Z(i>ph(ui) -f 2{TZ'^^Xir^^rn-l + K^_,) log(ZQ^) 



+2<°Li..® log(l +>® ) - 2 log - 2 ^ log 



+ -^3 



j = l 3 

where *p.i,.. stands for principal value integration, 7im is 
the dressing phase in the physical kinematics while i$ph 
is the same as i$ but evaluated in the physical region 
(see appendix). We used l/y^(wj ± i/2) = (following 
from (4)) and denoted Z[^^ = i^ + Yajl'^ ~ l/m^){u^ - 
y2'j5„,,2 isolate the poles in at u = Uj and to ensure 
decreasing asymptotics at large u which is of course useful 
for the numerics. Finally, in contrast to ^^^'^(ui), the 
term [ui] is evaluated in mirror kinematics. 



use the large L, ABA solutions of the y-system [14]. At 
each step of iterations we also update the position of the 
Bethe roots by solving the exact Bethe equation of the 
previous section. It is important to note that the r.h.s. 
of (5) remains, within our precision, purely imaginary in 
the process of iterations. 

We should also truncate the infinite set of y-functions 
to some finite number. We explicitly iterate the first 25 
y^ 's and 25 YqJs at each step. Then we also extrap- 
olate them to obtain extra 40 y^ 's and 40 Yq 's with 
higher numbers to replace the infinite sums in (4) on 
the next step of iteration. For Y*^ we can truncate the 
sums much earlier (typically Y^^/Y^^ <^ 0.1): the first 
5 of them are largely enough for our precision. Finally, 
the integrals along the real axis are computed along the 
stretch (—X, X) with X being a large cut-off. With these 
approximations, our absolute precision for the energy is 
around ±0.001. 

Wc solved the integral equations for a wide range of 
couplings < A < 700 ^ 1 stretching from the pertur- 
bative region up to this value, which is already a deep 
strong coupling regime [43]. We found no sign of any 
singularity and the curve looks perfectly smooth. By 
this reason we believe that any new singularities, such 
as those related to the Liischer /i-terms, will not appear. 
This seems to be the case perturbatively [37] and our 
numerics seems to indicate that the integral form of Y- 
system we are solving describes exactly the full spectrum 
of planar N=4 SYM theory in SL{2) sector. Although we 
cannot discard a possibility that some new singularities 
could collide with the integration contours for even larger 
values of the 't Hooft coupling (such extra singularities 
could be easily incorporated into our code), our numeri- 
cal results suggest that this possibility is very unlikely. 

With a precision of 0.001 we can approximate the 
Konishi dimension in the range we considered by 

4 A^rT~r 252. 93/1 -1-384. 74/1^-1-674. 13h^ + 128.17h-|-4 „,-u„^„ /, _ 
V9 + ^ 35. 67/1^-1-51. 43/1^-1-99.71/1^-1-29. 29/t+l wneit, IL — 

I \J -t- 1 and A = IGtt^ji^ (this function is just a short- 
hand for the data in Table I). 



CONCLUSIONS 



NUMERICS AND ITS INTERPRETATION 

We solve the integral equations (4) iteratively for the 
Konishi state. Namely, we find the y-functions at a step 
n by plugging the y-functions of step n — 1 into the r.h.s. 
of these equations. As the first step of the iterations we 



We presented here a numerical method for solving the 
Y-system for the Konishi dimension in planar N=4 SYM. 
It opens the way to the systematic study of the spectrum 
of many interesting states at any values of the coupling. 
We also hope to simplify the Y-system in the future using 
the underlying Hirota intcgrable dynamics and reducing 
it to a finite system of integral equations. 
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APPENDIX 



(w). For Konishi L = 2 and it has M = 2 derivatives. 



We use: r^-\v) ^ -^£log^, 



2n/ TT 1 

= 2 1^2 a,nd 



2i 



Q -, / [+n] _ [-m] 

(u, «) = ^ log , , - ^T'^") («) 



2i 



n— 1,771— 1 ' 



2TTi 

= T:oT:^oK2{u - V), Ktrn = K°F>K^{u ~ v) , 

where the fusion operation J"" o ^ = ^ A{u + ik). 

k— — 

Finally, we also use a nice integral representation [15] of 
the kernel %i,m given by 



CO „ 

-2 ^ y [S^A"^ («, - + ) Bt'J {w-^lv)+ C.C. 



(6) 



where sii? = j;'^ o j;^; o For the exact Bethe 

equations we should use this kernel in the mixed repre- 
sentation, 

CJO 

a=l 

Finally the source term = J^" o with 



(f>(w) = - log 

i 



(7) 



the 



BES [13] dressing phase S{u) ~ 
should be taken in the mixed, 



X 



where 

mirror-physical representation in the integral equations 
and in the physical-physical representation for $ph 
appearing in the exact Bethe equations. We use the 
mirror-physical integral representation [15] 



log S = log 



(8) 



iQ) * log 



R^+\u-iQ) 
R(-){u-iO) 



+ c.c. 



where (u) = ^ log p|i[_^_™j while for the physical- 
physical representation we can use the DHM integral 
representation [38]. Finally /Cph,mr(u, ii) represents a 
kernel where we use the physical (mirror) branches for u 
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Table I: Numerical values of the Konishi anomalous dimen- 
sion. The numerical error should be within ±0.001 
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6.3342 
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9.1149 


2.05 


10.5265 


0.50 


5.7120 






1.55 


9.2519 







[1] 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 



J. M. Maldacena, Adv. Theor. Math. Phys. 2 (1998) 
231 [Int. J. Theor. Phys. 38 (1999) 1113]; S. S. Gub- 
ser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 
428 (1998) 105; E. Witten, Adv. Theor. Math. Phys. 2 
(1998) 253 

J. A. Minahan and K. Zarembo, JHEP 0303 (2003) 013 
I. Bena, J. Polchinski and R. Roiban, Phys. Rev. D 69 
(2004) 046002 

N. Beisert, C. Kristjansen and M. Staudacher, Nucl. 
Phys. B 664 (2003) 131 

V. A. Kazakov, A. Marshakov, J. A. Minahan and 

K. Zarembo, JHEP 0405, 024 (2004) 

N. Beisert, V. A. Kazakov, K. Sakai and K. Zarembo, 

Commun. Math. Phys. 263, 659 (2006) 

R. A. Janik, Phys. Rev. D 73, 086006 (2006) 

M. Staudacher, JHEP 0505, 054 (2005) 

N. Beisert, Adv. Theor. Math. Phys. 12, 945 (2008) 

G. Arutyunov, S. Frolov and M. Staudacher, 10 

(2004) 016. 

N. Beisert, R. Hernandez and E. Lopez, JHEP 0611 
(2006) 070 

M. Bianchi, S. Kovacs, G. Rossi and Y. S. Stanev, JHEP 
0105 (2001) 042.; B. Eden, C. Jarczak, E. Sokatchev and 



5 



Y. S. Stanev, Nucl. Phys. B 722, 119 (2005). 
[13] N. Beisert, B. Eden and M. Staudacher, J. Stat. Mech. 

0701 (2007) P021 
[14] N. Gromov, V. Kazakov and P. Vieira, arXiv:0901.3753. 
[15] N. Gromov, V. Kazakov, A. Kozak and P. Vieira, 

arXiv:0902.4458 [hep-th]. 
[16] D. Bombardelli, D. Fioravanti and R. Tateo, 

arXiv:0902.3930 [hep-th]. 
[17] G. Arutyunovand S. Frolov, arXiv:0903.0141 [hep-th]. 
[18] AL B. Zamolodchikov, Phys. Lett. B 253 (1991) 391. 
[19] N. Dorey, J. Phys. A 39, 13119 (2006) 
[20] M. Takahashi, Cambridge University Press, 1999. 
[21] F.H.L. Essler, H.Frahm, F.Gohmann, A. Kliimper and 

V. Korepin, Cambridge University Press, 2005. 
[22] V.Bazhanov, S.Lukyanov and A. Zamolodchikov, Nucl. 

Phys. B 489, 487 (1997) 
[23] P. Dorey and R. Tateo, Nucl. Phys. B 482, 639 (1996) 
[24] D. Fioravanti, A. Mariottini, E. Quattrini and F. Ra- 

vanini, Phys. Lett. B 390 (1997) 243 
[25] G. Arutyunov and S. Frolov, JHEP 0903 (2009) 152 
[26] J. Ambjorn, R. A. Janik and C. Kristjansen, Nucl. Phys. 

B736 (2006) 288-301 
[27] G. Arutyunov and S. Frolov, JHEP 0712 (2007) 024 
[28] S. Frolov and R. Suzuki, arXiv:0906.0499 [hep-th]. 
[29] Z. Bajnok and R. A. Janik, Nucl. Phys. B 807, 625 (2009) 
[30] Z. Bajnok, R. A. Janik and T. Lukowski arXiv:0811.4448 
[31] F. Fiamberti, A. Santambrogio, C. Sieg and D. Zanon, 

Nucl. Phys. B 805 (2008) 231; V. N. Velizhanin, 

arXiv:0811.0607 [hep-th] 
[32] N. Beisert, T. McLoughlin and R. Roiban, Phys. Rev. D 

76 (2007) 046002 [arXiv:0705.0321 [hep-th]]. 
[33] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Nucl. 



Phys. B 636 (2002) 99 
[34] R.Roiban, A.Tseytlin, arXiv:0906.4294, (see 

A. Tseythn's talk at Shifmania 2009). 
[35] A. Tirziu and A. A. Tseytlin, Phys. Rev. D 78 (2008) 

066002 

[36] V. Kazakov, A. Sorin and A. Zabrodin, Nucl. Phys. B 
790, 345 (2008) 

[37] Z. Bajnok, A. Hegedus, R. Janik, T. Lukowski, 
arXiv:0906.4062 [hep-th]. 

[38] N. Dorey, D. M. Hofman and J. M. Maldacena, Phys. 
Rev. D 76 (2007) 025011 

[39] Here Z is one of the complex scalars and D is a covariant 
derivative in a light cone direction. 

[40] L is the number of Z fields in an operator in SL{2) sector. 

[41] meaningful until the convergency radius g < gc = 1/4; 
in this region our numerical data are indistinguishable 
from the BAE for our accuracy. At week coupling, an 
interesting numerical prediction could be the order g^". 
However our numerical error ±0.001 is larger than (t/c)'^" 
and therefore does not allow for such predictions. 

[42] Assuming the leading coefficient to be 2X^^'^ one gets 
1.001 instead of 0.994 for the subleading term, even closer 
to 1 - the value predicted here. 

[43] The expected world-sheet perturbation theory expansion 
at strong coupling is A^^*(ai -I- a^j^fX -\- . . . ). The first 
two coefficients - are of order 1 and we assume that all 
ai ~ 1. In our numerics 1/-v/A ~ 0.04. The appropriate 
extrapolation procedure can easily increase the precision 
by one order and we expect the absolute numerical error 
for our 02 coefficient to be within 0.005. 



